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We study color superconductivity with Nf — 1,2, and 3 massless flavors of quarks. We present a 



general formalism to derive and solve the gap equations for condensation in the even-parity channel. 
This formalism shows that the leading-order contribution to the gap equation is unique for all color 
OO ' superconductors studied here, and that differences arise solely at the subleading order. We discuss 

a simple method to compute subleading contributions from the integration over gluon momenta 
in the gap equation. Subleading contributions enter the prefactor of the color-superconducting 
gap parameter. In the case of color-flavor and color-spin locking we identify further corrections to 
this prefactor arising from the two-gap structure of the quasiparticle excitations. Computing the 
transition temperature, Tc, where the color-superconducting condensate melts, we find that these 
contributions lead to deviations from the BCS behavior Tc — 0.57 00, where 0o is the magnitude of 
. the zero-temperature gap at the Fermi surface. 

o ■ 

O ■ I- INTRODUCTION AND CONCLUSIONS 



> 

o 
in 
o 



Cold and dense quark matter is a color superconductor [1,2]. At asymptotically large quark density or, equivalently, 
quark chemical potential /i, asymptotic freedom [3] implies that the strong coupling constant g becomes small. In this 
case one can reliably compute the color-superconducting gap parameter to leading and subleading order in g from 
(-H J a gap equation derived within the framework of QCD [1,4,5]. For instance, in a color superconductor with Nf = 2 
massless flavors of quarks (commonly called the "2SC" phase), the value of the gap at the Fermi surface and at zero 
temperature is 
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.fc^266[,A.exp(-^j , (1) 
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where 

The term in the exponent of Eq. (1) was first computed by Son [6]. It arises from the exchange of almost static 
magnetic gluons. The factor b in front of the exponential originates from the exchange of static electric and non-static 
magnetic gluons [4,5]. The prefactor fej, is due to the quark self-energy [7,8]. 

In color superconductors with Nf = 1 and 3 flavors, various other prefactors may arise [9,10], but the exponential 
exp[— 7r/(2^)] remains the same. As will be demonstrated in this paper, this is not an accident, but due to the fact 
that the leading-order contribution to the QCD gap equation does not depend on the detailed color, flavor, and 
Dirac structure of the color-superconducting order parameter. This structure only enters at subleading order, and we 
provide a simple method to extract these subleading contributions. 

Let us briefly recall what the terms "leading," "subleading," and "sub-subleading order" mean in the context of 
the QCD gap equation [8]. Due to the non-analytic dependence of 4>o on the strong coupling constant g one cannot 
apply the naive perturbative counting scheme in powers of g in order to identify contributions of different order. In 
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the QCD gap equation there are also logarithms of the form ln(/i/(/)o), which are ^ \/g due to Eq. (1) and thus 
may cancel simple powers of g. A detailed discussion of the resulting, modified power-counting scheme was given in 
the introduction of Rcf. [8] and need not be repeated here. In short, leading-order contributions in the QCD gap 
equation are due to the exchange of almost static magnetic gluons and are proportional to g^^o ln^(A*/^o) ~ '/'o- 
They determine the argument of the exponential in Eq. (1). Subleading-order contributions are due to the exchange 
of static electric and non-static magnetic gluons and are ^ g"^ (pQ ln(/i/(/>o) ^ g 4>o. They determine the prefactor of the 
exponential in Eq. (1). Finally, sub-sublcading contributions arise from a variety of sources and, at present, cannot 
be systematically calculated. They are proportional to .g^ c/io and constitute 0{g) corrections to the prefactor in Eq. 
(1). It was argued that also gauge-dependent terms enter at this order [11]. This is, of course, an artefact of the 
mean- field approximation which was used to derive the QCD gap equation [12]. On the quasiparticle mass shell, the 
true gap parameter is in principle a physical observable and thus cannot be gauge dependent. 
In color superconductors, the mass shell of a quasiparticle is determined by its excitation energy 

ekA<l>) = [{k - l^f + K \H^k,r, k)f] , (3) 

where k = |k| is the modulus of the 3-momentum of the quasiparticle, and (f){€k,r,k) is the gap function on the 
quasiparticle mass shell. The index r labels possible excitation branches in the superconductor, which differ by the 
value of the constant A,-. At the Fermi surface, k = fi, the excitation of a quasiparticle - quasiparticle-hole pair costs 
an energy 2e^,r(^^) = 2\/A^^o- The true energy gap is therefore \/X^(f)o- 

At first sight the introduction of the constant appears somewhat awkward. The advantage is that it allows 
to generalize Eq. (3) to different color-superconducting systems. For example, in a two-flavor color superconductor, 
quarks of two colors form Cooper pairs with total spin zero, while the third color remains unpaired [1] . Consequently, 
there are two difi^erent excitation energies, ek^i and ek^2- Four quasiparticle excitations have Ai = 1, with gap (j)o, 
while two have A2 = 0, corresponding to the unpaired quarks. These are so-called "ungapped" excitations. At the 
Fermi surface, it costs no energy to excite them. In a three-flavor color superconductor, with color-flavor locking 
(CFL) [13], all nine quark colors and flavors form Cooper pairs, but there are still two distinct branches of fermionic 
excitations. The first, with Ai = 4, occurs with degeneracy one, while the other, with A2 = 1, has degeneracy eight. 
The gap corresponding to the first excitation has magnitude 2(f)Q, while for the other eight the size of the gap is (j)Q. 
A similar two-gap structure also appears in the color-spin locked (CSL) phase of a one-fiavor color superconductor 
[10]. However, here the first excitation, with Ai = 4, has a four-fold degeneracy, while the second, with A2 = 1, has 
an eight-fold degeneracy. 

In this paper we aim to clarify the similarities and differences between various color-superconducting systems. To 
this end we systematically study six different cases. The first two cases are spin-zero color superconductors with (i) 
two flavors of massless quarks and (ii) three flavors of massless quarks in the CFL phase. The other four cases deal 
with one massless quark flavor. In this case, the condensate has spin one [14]. Similar to Helium-3 [15], this allows 
for a multitude of different phases, distinguished by the symmetries of the order parameter [10]. We only focus on 
(iii) the CSL phase with longitudinal and transverse gaps, (iv) the CSL phase with longitudinal gap only, (v) the 
CSL phase with transverse gap only, and (vi) the polar phase. In this context, "longitudinal" and "transverse" refers 
to pairing of quarks with the same or different chiralities, respectively [5]. The cases (iv) and (v) can be considered 
separately because, as we shall show below, longitudinal and transverse gaps do not induce each other. In all six cases 
we only consider condensation in even-parity channels, because these are favored by effects which explicitly break the 
U{1)a symmetry of the QCD Lagrangian. 

We show that in all six cases the gap equation has the general form 

<^(efe,., k) = f f d{q - M) E ^(^..«) t^^l^ iw) \ f ^^^1 • (4) 

The sum over s runs over all distinct branches of fermionic excitations with energy Cq^s in the color superconductor. 
For the systems considered here, there are only two such branches, such that s = 1 or 2. The coefficients are 
positive numbers, obeying the constraint 

2 

E«« = l- (5) 

In the first four columns of Table I we display the values of A^ and for the six cases studied here. There are 

ungapped excitations (A ^ = 0) in the cases (i) , (iv) , (v) , and (vi) , while all excitations are gapped in the cases (ii) and 
(iii). For ungapped excitations, the corresponding Os vanishes, and thus these do not appear in the gap equation. This 
is natural, because ungapped excitations should not affect the value of the color-superconducting gap. The constants 
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are identical for the cases (ii) and (iii), Ai = 4 and A2 = 1. The coefficients also assume the same values, 1/3 and 
2/3, but the association of these values with the corresponding constants Ag is reversed in case (iii) as compared to 
case (ii). 

The occurrence of the wave function renormalization factor Z{eq^s) in a gap equation of the type (4) was first 
discussed in Ref . [8] . The constant b in Eq. (4) is defined as 

6 = 6 exp(— d) , (6) 

with b from Eq. (2), and d a constant of order one. The constant d originates from subleading contributions to the 
gap equation. For spin-zero condensates, d = 0, due to an accidental cancellation of some of the subleading terms 
arising from static electric and non-static magnetic gluon exchange. In the spin-one cases, this cancellation does not 
occur and, consequently, d ^ 0. 

In this paper we present a simple method to extract the value of the constant d without actually solving a gap 
equation. This method utilizes the fact that, to subleading order, the integration over gluon momenta in the QCD 
gap equation can be written as a sum of a few integrals multiplied by constants. Only these constants depend on the 
detailed color, flavor, and Dirac structure of the order parameter. The integrals are generic for all cases studied here 
and have to be computed only once. The precise numerical values for d are listed in the fifth column of Table I. 

The fact that we can write the gap equation in all six cases in the form (4) is nontrivial. It means that the leading 
contribution to the gap equation is unique. If it were not, then the prefactor of the gap integral would be different 
for each case. In other words, the contribution of almost static magnetic gluons to the gap equation is universal in 
the sense that it is independent of the detailed color, fiavor, and Dirac structure of the color-superconducting order 
parameter. Differences between the six cases studied here occur at subleading order. Only at this order the specific 
structure of the order parameter is important and leads to different values for the constant d in Eq. (6). 

We solve the gap equation (4) at zero temperature and compute the value of the gap function at the Fermi surface, 
(^0- In all cases studied here, we can write the result in the form 

<^o = 2bb',n exp (^-^) (Ar A^)-'/' • (7) 

(Remember that 0*^ = 1.) Prom this equation and Eq. (1) one immediately determines ^0 in units of the gap in the 
2SC phase, 

-^=exp(-d) (X^XT)-'^' . (8) 

This ratio is given in the sixth column of Table I. For spin-one color superconductors, d is positive, and the exponential 
factor leads to a tremendous suppression of the gap by factors e^^'^ ~ 10~^ to ~ 2.5 x lO^'^ relative to the spin- 
zero gap [7,10]. In contrast to the value of d, the additional factor (A^^Ag^)"^^^ in Eq. (8) cannot be simply read off 
from the subleading contributions in the gap equation, but only follows from the explicit solution. It is different from 
1 for color superconductors with two distinct branches of gapped quasiparticle excitations. In this case, this factor 
further reduces ^0 as compared to the 2SC case. 

1 /2 

The factor {X^'X^^y ' is also different from 1 in case (v) where Ai = 2. However, according to Eq. (3) the "true" 
gap is \/Xi4'o = V^4'o, and not 0o- Therefore, the ratio of the true gap to the gap in the 2SC case is just exp(— d). 
In order to indicate this, in Table I we put the factor 2~^/^ arising from Eq. (8) in parentheses. 
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2 





i 





l)/2 


(2 " 




(vi) 


polar 


1 





1 





1(3-1- cos" 1?) 




1 



TABLE I. The constants As, as, d, the ratio 0o/(io" % and the ratio Tc/(j)o normalized to its BCS value. In case (vi), ■& is the 
angle between the direction of the spin-one condensate and the 3- momentum of the quarks in the Cooper pair, see Sec. 11 G. 
In case (v), the factors in parentheses do not occur if 0o is replaced by the true gap \/\i4>o- 
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Finally, we discuss the transition temperature Tc, where the color-superconducting condensate melts. We find 

^ = £l (A^i Xl^fl^ ~ 0.57 Xl^fl'' , (9) 

where 7 ~ 0.577 is the Euler-Mascheroni constant. In the cases (i), (iv), and (vi), where there is only one gapped 
quasiparticle excitation, (A°^ Aj^)^/^ — 1, and we recover the relation Tc/(/)o — 0.57 well known from BCS theory 
[16]. Its validity for QCD with Nj = 2 flavors of massless quarks was first demonstrated in Refs. [5,8]. In case (v), 
(Ai^ Aa')^/^ = y/2, but this factor is absent if we rescale in Eq. (9) by to obtain the true gap. Therefore, also 
in this case the BCS relation between the zero-temperature gap and the critical temperature is valid. In the cases 
(ii) and (iii) there are two distiiic;t gapped quasiparticle excitations, and consequently two gaps, •v/Ai(/>o = 2(/>o and 
vl^^o = (t>o- The BCS relation Tc/(po = c'^/tt is violated by the additional factor (AJ^ A^')^^ > i 

In order to elucidate the deviations from the BCS relation, in the last column of Table I we present our results for 
Tc in units of the critical temperature expected from BCS theory. Apparently, the two-gap structure in the cases (ii) 
and (iii) is responsible for the observed deviations. It would be interesting to observe similar behavior in other weak- 
coupling superconductors with more than one gapped excitation branch. Note that, although Tc/<j)o is different than 
in BCS theory, the absolute values of do not change. If the energy scale is set by ^o^*"^, then Tf./(j)'^'~' = exp(— d), 
because the factor (A^^ Aa')"^/^ in Eq. (8) simply cancels the factor (A"^ X'^'Y^'^ in Eq. (9). 

The remainder of this paper is organized as follows. In Section II we show that the gap equations for all six 
cases considered in this paper is of the form (4). We explain the origin of the constants A^ as eigenvalues of an 
operator constructed from the color-superconducting gap matrix. We also present a simple method to compute 
subleading corrections to the gap arising from the integration over gluon momenta in the gap equation, leading to 
the suppression factor exp(— d) in Eq. (6). In Section III we solve the gap equation (4) at zero temperature and 
explain the occurrence of the additional factor (A^^ Aj^)"^/^ in Eq. (7). Finally, in Section IV we compute the critical 
temperatures T^.. 

Our convention for the metric tensor is g^'^ = diag{l, —1, —1, —1}. Our units are h = c = ks = Four- vectors are 
denoted by capital letters, K = = (fco, k), and k = |k|, while k = k/fc. We work in the imaginary-time formalism, 
i.e., T/V^i^ = TJ2n I rf^k/(27r)^, where n labels the Matsubara frequencies a;„ = iko. For bosons, a;„ = 2n7rT, for 
fermions, u>n = (2n -|- l)7rT. 



II. GAP EQUATIONS 
A. General derivation 

In fermionic systems at non-zero density, it is advantageous to treat fermions and charge- conjugate fermions as 
independent degrees of freedom and to work in the so-called Nambu-Gorkov basis. In this basis, the full inverse 
fermion propagator is defined as 

^-1 ^ f "Sill 5*12! \ ^ f S^i^ + Ell S12 I ^ ^^Q^ 

where S^u is the propagator for free fermions, iS''^22 the propagator for free charge- conjugate fermions. In momentum 
space and for massless quarks, 

S\iiK) = irX^ + /x7o)"' , S°22{K) = i^K^ - ^70)-' , (11) 

where ■y^ are the Dirac matrices. The 11 component of the self-energy, En, is the standard one- loop self-energy for 
fermions; similarly, S22 is the self-energy for charge-conjugate fermions. In Ref. [8] it was shown that, in order to 
solve the gap equation to subleading order, it is permissible to approximate these self-energies by 

S(if) = En(i^) = E22(if) 70 9^ koln^ , (12) 

where = {3TT/4)ni^; the zero-tcmpcrature gluon mass parameter (squared) is — N f g'^ /j? / . The 21 
component of the self-energy, E21, which was denoted in [5], is the gap matrix in a superconductor, while 
E12 = 7oE^i7o. 

Inverting Eq. (10) one obtains the full fermion propagator S. The 11 component. 
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'S'll — (S''^22 +^22) {S^n +^ll)('S'°22 +S22)— C 

is the full quasiparticle propagator, where we defined 

C = El2(S'°22 + 5]22)~^5]2l(S'°22 +^22) • 

The 21 component is the so-called "anomalous" propagator. It is given by 

521 = — ('S''^22 + 5^22)"^ ^21 Sii . 

In all cases considered here, the gap matrix can be written as 

S2i(i^) = ^</.^(if)7WkA£, 

e=± 



(13) 



(14) 



(15) 



(16) 



where cj)'^ {K) is the gap function, A^k is a matrix defined by the symmetries of the color-superconducting condensate, 
and = (1 + 6707 • k)/2, e = ±, are projectors onto states of positive or negative energy. In general, A^k is a matrix 
in color, flavor, and Dirac space, and is constructed such that 



[A4k,A£] = 0. 

With the gap matrix (16), the operator C{K) assumes the form 

C(/^) = ^|0^(if)pLkA 



k ' 



(17) 



(18) 



where 



= 70 M\.Mk 70 • 



(19) 



Note that also [Lk, A^] = 0. Since Lk is hermitian, it has real eigenvalues and can be expanded in terms of a complete 
set of orthogonal projectors V^, 



k : 



(20) 



where \r are the eigenvalues of Lk- Our choice of the symbol is judicial: it will turn out that they are identical with 
the constants \r appearing in the quasiparticle excitation energy (3) and which are listed in Table I. In Appendix A 
we determine the eigenvalues of Lk and their degeneracy for the six color-superconducting systems studied here. 

In all cases considered in this paper, there are only two distinct eigenvalues, so that one can easily express the two 
corresponding projectors in terms of Lk, 



■>1,2 



-^2.1 



Al,2 — ^2,1 



(21) 



Obviously, these projectors also commute with the energy projectors, [P^' , A^] = 0. 

The next step is to compute the full quasiparticle propagator ^n. The inversion of the term in brackets in Eq. (13) 
is particularly simple, because the four projectors "Pk^A^ are orthogonal and form a complete set in color, flavor, and 
Dirac space. With Eqs. (12), (13), (18), and (20) we obtain 



S22 



1 



k-"^k 



(22) 



where 



Z{ko) 



1 2, 



(23) 



is the wave function renormalization factor introduced in Ref. [17] and 
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21 1/2 



(24) 



arc the excitation energies for quasiparticles, e = +, see Eq. (3), or quasi- antiparticles, e = — . 

In order to compute the anomalous propagator 521, we insert Eq. (22) into Eq. (15) and employ Eq. (16). The 
result is 



In the mean-field approximation [12], S21 obeys the gap equation [8] 



(25) 



(26) 



where are the Gell-Mann matrices (times a factor 1/2) and A"^ is the gluon propagator. 

To derive the gap equation for the gap function cff^K), we insert Eq. (25) into Eq. (26), multiply both sides from 
the right with and trace over color, flavor, and Dirac space. To subleading order in the gap equation, it is 

permissible to use the gluon propagator in the Hard-Dense-Loop (HDL) approximation [18], where it is diagonal in 
adjoint color space, A^^ = Sab^'^"- We obtain 



A''''(X-g)r--«(k,q) 



(27) 



where 



^;^^(k,q) 



Tr 



Tj 70 70 Aq 7, Ta A^ 



Tr 



(28) 



The form (27) of the gap equation holds for all cases considered in this paper. What is different in each case is the 
structure of the term T^^ ■■'(k, q). Our computation will be done in pure Coulomb gauge, where 



A""(P) = A^(P) , A"'(P) = , A'^(P) = [5^^ - fp>) At{P) 



(29) 



with the longitudinal and transverse propagators A^^t and P = K — Q. Consequently, we only need the 00-component, 
Tqq '*(k, q), and the transverse projection of the ij-components. 



rr ■^(k,q)^-(<5^-' -pV)^r ''(k,q) , 



(30) 



of the tensor (28). (The extra minus sign is included for the sake of notational convenience.) It will turn out that in 
all cases studied here the quantities TQ^^^^k, q) are related in the following way: 



7;,^r^(k.q) _7;-'2(k,q) 



r-'^(k,q) r-'^(k,q) 



const. . 



(31) 



The right-hand side of Eq. (28) depends on fc, q, and k • q. The latter can be replaced by the square of the gluon 
3-momentum via k • q = (A;^ -|- 9^ — p'^)/{2kq). Thus, the relevant components can be written in terms of a power 

series in , 



To^o^''«(k,q)=a« ^ vLiee',k,q) 



T-''«(k,q)=a, E Vlm{ee',k,q) f 



m——l 



2 - 

2 \ ™ 



kq 



(32a) 
(32b) 
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Here, the coefficients rl-^mk^^ ^ longer depend on s on account of Eq. (31). The overaU normahzation on the 

right-hand side of Eq. (32) is still free, and we choose it such that Eq. (5) is fulfilled. This uniquely determines the 
values of the dimensionless coefficients rl^ljyee! , k, q). 

We now perform the Matsubara sum in Eq. (27), which does not depend on the detailed structure of the tensor 
q-ee ,s(ij;^q^_ This calculation is similar to that of Ref. [5]. The difference is the appearance of the wave function 
renormalization factor Z{qo) [8]. To subleading order, this amounts to an extra factor Zle^ g) in the gap equation. 
Since there are two different excitation energies Cq^i and eq 2 on the right-hand side of the gap equation, we can put 
the gap function on the left-hand side on either one of the two possible quasiparticle mass shells ko = ek,i or ko = ek,2- 
One then obtains 
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dqqY,asZ{ei,) 



tanh 



+ 



^Q{p-M) + Q{M-p) 



+ ■ 



V2n 



(33) 



The first term in braces arises from static electric gluons, while the two terms in brackets originate from non-static 
and almost static magnetic gluons, respectively. Various other terms which yield sub-STibleading contributions to the 
gap equation [5] have been omitted. In deriving Eq. (33) we assumed that the gap function does not depend on the 
direction of k. This is true in all cases considered here, except for the polar phase, where we neglect this dependence, 
cf. Sec. II G. 

Although the coefficients 772^ depend on k and q, to subleading order in the gap equation we may approximate 
k q ^. This can be easily proven by power counting. To this end, it is sufficient to take k = fi, and write q = + 
where ^ = g — A*. In weak coupling, the gap function is sharply peaked around the Fermi surface, and thus the range 
of integration in the gap equation can be restricted to a small region of size 2 5 around the Fermi surface. All that 
is necessary is that 5 is parametrically much larger than i/io, but still much smaller than /i, 4>o S 1.1 [5]. It turns 
out that 6 ~ nig is a convenient choice. Since the integral over ^ is symmetric around ^ = 0, terms proportional to 
odd powers of ^ vanish by symmetry. Thus, corrections to the leading-order terms are at most ~ (^/yu)^. As long as 



S is parametrically of the order of rrig, ^ < nig, and these corrections are 



suppressed by two powers of the 



coupling constant. Even for the leading terms in the gap equation the correction due to terms ~ (^/a*)^ is then only 
of sub-subleading order and thus negligible. 

Since the coefficients 772m dimensionless, with the approximation k c^l q ji they become pure numbers which, 
as we shall see in the following, are directly related to the constant d discussed in the introduction and listed in Table 
I. In all cases considered here, T^jm = ^^^^ m > 3, and the series in Eq. (32) terminate after the first few terms. 
Moreover, rj'^_2 always vanishes and, to subleading order, also rf_2 = 0. For the remaining m, the p integral in Eq. 
(33) can be performed exactly. The details of this calculation are deferred to Appendix B. We obtain 



''^^lr,k) = ^ ^ dqY,asZ{el,) 



r'(<„9) 



tanh 




Vol 111 



4„2 42 



(34) 



Note that the contribution from almost static magnetic gluons only appears in the term proportional to tjq, while 
non-static magnetic and static electric gluons contribute to all other terms. 

The antiparticle contribution (e' = — ) does not have a BCS logarithm, since e~g ~ g -|- /i. For the same reason, 
for antiparticles the logarithm from almost static magnetic gluons is also only of order 1, and furthermore there is no 
large logarithm from the p integrals. Therefore, the antiparticles contribute at most to sub-subleading order to the 
gap equation and can be neglected. In the following, we may thus set e = e' = -|- and omit this superscript for the 
sake of simplicity. Then, the gap equation for the quasiparticle gap function reads 



Hek,r, k) = f 1^ d{q - M) E ^(^«.«) ta^li i^) i ^0 ^ 



Ie2 -e2 



(35) 



k.r I 



where 



Ml/ 



exp(— 2d) 



(36) 
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with 

d=-^ H+vl + 2ivi + vl)] ■ (37) 

In all cases considered in this paper, rj^ = rj^, so that b assumes the value quoted in Eq. (6). The expression (37) is a 
general formula to compute the constant d from the coefficients r/j^^- We also find that, for all cases considered here, 
ril = 2/3. This is the uniqueness of the leading-order contribution to the gap equation mentioned before. With this 
value of ?7q, the gap equation has the general form (4). 

In the following, we shall discuss spin-zero color superconductors in the 2SC and CFL phases, as well as spin-one 
color superconductors in the CSL phase with both longitudinal and transverse gaps, the CSL phase with longitudinal 
gap only and with transverse gap only, and the polar phase. Each case is uniquely characterized by the matrix A^k 
which is given by the symmetries of the color-superconducting condensate. This matrix determines the eigenvalues 
and the projectors P^. Evaluating the traces in Eq. (28) and comparing with Eq. (32), one reads off the coefficients 
ri2^, as well as the constants o^. This completely specifies the gap equation in each case. 



B. The 2SC phase 



For Nf = 2, the spin-zero condensate is a singlet in flavor and an antitriplet in color space [1]. The (antisym- 
metric) singlet structure in flavor space can be represented by the second Pauli matrix (r2)/g = ie/g, f,g = 1,2. 
The (antisymmetric) antitriplet structure in color space restricts the gap matrix to be a linear combination of the 
antisymmetric Gell-Mann matrices A2, A5, and Ay. These form an 5*0(3) subgroup of SU{3)c, so that we can also 
choose the generators {Ji)jk = —i^jki ii3, k = 1,2, 3, of 5*0(3). The gap matrix is thus a scalar in flavor space and 
a 3-vector in color space. Upon condensation, this vector points in an arbitrary, but fixed, direction which breaks 
SU{2)c to SU{2)c. For the sake of convenience, we align this vector with J3. Thus, the matrix A^k reads 



A^k = h T2 75 



(38) 



where 75 takes into account that we restrict our discussion to the even-parity channel. This matrix obviously fulffils 
the condition (17). From Eq. (19) we construct the matrix 



(39) 



In this case, Lk does not depend on k, and consists of a unit matrix in flavor space and a projector onto the first two 
colors in color space. In principle, it also consists of a unit matrix in Dirac space, which we disregard on account of 
the spin-zero nature of the condensate. 

The eigenvalues of Lk are (cf. Appendix A) 



Ai = 1 (4-fold) 



A2 = (2-fold) 



From Eq. (24) we conclude that there are four gapped and two ungapped excitations. 
The projectors Vk follow from Eq. (21), 



(40) 



(41) 



They have the property that J^P^ 
s = 1 we obtain 



^oo'''(k,q) 



Tr'''(k,q) 



J3 and Js'Pk = 0. Consequently, the tensor T^^ '^(k, q) vanishes trivially. For 

i (l + ee'k-q 
3 — ee' k • q 



(efc - e'qf 

p2 



(l + ee'k-q) 



(42a) 
(42b) 



We now match this result to the expansion in terms of p^, Eq. (32). Since T^^ '^(k, q) = and because of Eq. (5), we 
have 



ai = 1 



02 = . 



(43) 
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This uniquely fixes the coefBcients ??2m(^^'' ^' l)- subleading order we only require their values for e = e' = + and 

A: ~ g ~ /z, 

2 1 2 1 

??o=3 > 4 = --^ , ??i = , ^0 = 3 > ^2=g , »?1 = 0. (44) 

This result implies that the contributions from static electric and non-static magnetic gluons to the constant d defined 
in Eq. (37) cancel, and consequently d = Q. 



C. The CFL phase 

In the CFL phase, the spin-zero condensate is a flavor antitriplet locked with a color antitriplet [13], 

A^k = J-l75, (45) 

where J = (Ji, J2, J3) represents the antitriplet in color space, with {Ji)jk = —i^ijk as introduced above. The vector 
I represents the antitriplet in flavor space and is deflned analogously. Consequently, (J • = —S( Sj + S? Sj . This 
condensate breaks SU{3)c x SU{3)f to SU{3)c+f- 
From Eq. (19) we obtain the matrix 

{L^){f = [(J • I)^] J = S{ <5| + S., Sf^ . (46) 

As in the 2SC case, the operator Lk is independent of k, and we omitted its trivial Dirac structure. It can be expanded 
in terms of its eigenvalues and projectors as in Eq. (20), with (cf. Appendix A) 

Ai = 4 (1-fold) , A2 = 1 (8-fold) , (47) 

and 

(r^l! = 1^(51 , {ri){f = SijSf<>-ls{si, (48) 

where and correspond to the singlet and octet projector introduced in Ref. [19]. 

We now compute the relevant components of the tensor T^^ '*(k, q). Since the Dirac structure of A^k is the same 
as in the 2SC case, the dependence on k and q is identical to the one in Eq. (42). However, since the color-flavor 
structure is different, we obtain a non-trivial result both for s = 1 and s = 2, with different prefactors, 

■^(k,q) = iTo-''^(k,q) = i (l + ee'k-q) , (49a) 



T-'i(k,q) = lTr'^(k,q) = i 



3_ee'k-q-fc^ (l + ee'k-q) 



(49b) 



Obviously, the condition (31) is fulfilled. The coefficients r^jj^ remain the same as in Eq. (44), which again yields 
d = 0. However, the two-gap structure leads to the constants 

1 2 

ai = - , "2 = 3 • (50) 

In our treatment we have so far neglected the color-sextet, flavor-sextet gap which is induced by condensation in the 
color-antitriplet, flavor-antitriplet channel [20]. Such a color-flavor symmetric structure is generated in the anomalous 
propagator 521, even for the completely antisymmetric order parameter of Eq. (45). (This does not happen in the 
2SC case, where the color-flavor structure of 52i remains completely antisymmetric.) Consequently, it also appears 
on the right-hand side of the gap equation. The reason why it disappeared in our calculation is that we projected 
exclusively onto the antisymmetric color-flavor channel when we multiplied both sides of Eq. (26) with Ad]^ A^ and 
traced over color, flavor, and Dirac space. To be consistent, one should have started with an order parameter which 
includes both the symmetric and the antisymmetric color-flavor structures. In weak coupling, however, the symmetric 
gap is suppressed by an extra power of the strong coupling constant g [9] . This fact by itself is not sufficient to neglect 
the symmetric gap in the weak-coupling solution of the gap equation because, as explained in the introduction, this 
could still lead to a subleading correction which modifies the prefactor of the (antisymmetric) gap. One way to avoid 
this is a cancellation of the leading terms involving the symmetric gap in the gap equation for the antisymmetric gap. 
A more detailed investigation of this problem, however, is outside the scope of the present paper. 
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D. The CSL phase 



For condensation in the even-parity, spin-one channel the gap matrix reads (cf. Appendix C, see also Ref. [5]) 

E2i{K) = J24>''{K)-[k + 7A^)]K, (51) 

e=± 

where 7j_ (k) = 7 — 7 • k k and 7 = (7^ , 7^ , 7'^) ■ This most general form for the gap matrix in the spin-one case differs 
from the one in Ref. [10] by the appearance of 7j^ (k) instead of 7. Prom the discussion in Appendix C it is obvious 
that both forms are equivalent. 

The spin-one condensate is an SU{2) triplet, and thus the order parameter (j)'^{K) is a 3- vector. In the CSL phase, 

each spatial component of this vector is assigned a direction in color space, (x,y,z) {r,g,b). This breaks color 
SU{3)c and spatial SO{3) to an SO{3) subgroup of joint color and spatial rotations. The matrix reads 



A^k = J - [k + 7±(k)] . 



(52) 



This matrix fulfills the condition (17), due to the fact that Aj^ commutes with 7j^(k). Had we used 7 in Eq. (51), 
like in Ref. [10], this condition would have been violated and the general discussion presented in Sec. II A would not 
apply to the subsequent calculation. 
Prom Eqs. (19) and (52) we compute 



{L^r^g = 5'^ 5„0 + U 5^^ + 7i«-,(k) y 6^0 - 7i^.(k) 



(53) 



In contrast to the 2SC and CFL cases, this 12 x 12 matrix in color and Dirac space now explicitly depends on k. 
Nevertheless, its eigenvalues are pure numbers (cf. Appendix A), 



Ai = 4 (4-fold) 
The projectors follow from Eq. (21), 

1 



A2 = 1 (8-fold) . 



i-Plyip = - \ [k' Sa, + 7ia-y(k)] S,0 - ^{^^{k) 

Inserting these projectors and from Eq. (52) into Eq. (28) we obtain 

i Tg'\k, q) = T-''2(k, q) = ^ (1 + ee' k . q) [l + (1 + ee') k • q 
i Tr''HK q) = T-''2(k, q) = 1 {2 k . q (1 - ee' k . q) 



+ 



(efc- 



(l + ee'k-q) [l + (1 + ee') k • q] | 



(54) 

(55a) 
(55b) 

(56a) 
(56b) 
(56c) 



Comparing this to Eq. (49), the prefactor 1/2 now accompanies T^q^'^ instead of TQo,t^- Consequently, the constants 
ai and 02 exchange their roles compared to the CPL case, Eq. (50), 



ai = 



02 = 



(57) 



and, to subleading order, 

' 2 
3 ' 



V2 



7_ 
18 



1 

18 



5 
18 



r?* = . 



(58) 



According to Eq. (37), this yields d = 5. 

As in the CFL case, another condensate with a symmetric color structure is induced. This condensate belongs to 
the color-scxtct representation and, for Nj- = 1, necessarily carries spin zero. To identify this induced condensate, 
one has to explicitly analyze the color structure of 521. By analogy to the CFL case, we expect this condensate to be 
suppressed by a power of g compared to the primary spin-one, color-antitriplet condensate. Its contribution to the 
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gap equation could be of sub-subleading order, if there is a cancellation of the leading terms involving the spin-zero 
gap in the gap equation for the spin-one gap. A more detailed investigation, however, is beyond the scope of the 
present paper. 

In the following two subsections we study two special cases of the CSL color superconductor. The first is A^k ^ 
J • k, and the second is A^k ~ J • 'y±0^)- In the first case, the gap matrix commutes with the chirality projector 
'Pr,e = (1 i 75)72, and consequently only quarks of the same chirality form Cooper pairs. The ensuing gap was termed 
longitudinal gap in Ref. [5]. (It corresponds to the LL and RR gaps of Ref. [10].) In the second case, commuting the 
gap matrix with the chirality projector flips the sign of chirality, which indicates that the quarks in the Cooper pair 
have opposite chirality. This leads to the so-called transverse gap [5] (the LR and RL gaps of Ref. [10]). The reason 
why we study both cases separately is that a purely longitudinal gap matrix on the right-hand side does not induce a 
transverse gap on the left-hand side of the gap equation and vice versa. This will be explained in more detail below. 

E. The longitudinal CSL phase 

In the CSL phase with longitudinal gaps only, the matrix Aiu reads 

A4k = J • k . (59) 
The condition (17) is trivially fulfilled. Inserting Eq. (59) into Eq. (19), we obtain 

{Lk)% = {S'^ - y) Sap . (60) 

This matrix is a projector onto the subspace orthogonal to k. However, due to color-spin locking, the indices i,j run 
over fundamental colors and not over spatial dimensions, and thus, amusingly, this projection actually occurs in color 
space. Since Lk is a projector, we find the eigenvalues (cf. Appendix A) 

Ai = 1 (8-fold) , A2 = (4-fold) . (61) 
The projectors P^'^ follow from Eq. (21), 

^k = ik , 7'^ = 1-Lk, (62) 

similar to the 2SC case, cf. Eq. (41). The peculiar feature of Eq. (62) is that the projector belongs to the 
eigenvalue corresponding to quasiparticle excitations with a longitudinal gap, but it actually projects onto the subspace 
orthogonal to k. This is, however, not a contradiction, since the projection occurs in color space, while the gap is 
longitudinal (parallel to k) in real space. 

The similarity to the 2SC case carries over to the quantities T^Q^^{]s.,q). For s = 2, these quantities again vanish 
because J • kP^ = 0. For s = 1, we obtain 



(k,q) = ik q (l + ee'k-q) , (63a) 



T-''^(k,q) = ik.q [s-ee'k-q-^^^^^ (l + ee'k-q) 



(63b) 



which only diSev by an overall factor k • q from those of Eq. (42). While the constants a^. are the same as in the 2SC 
case, see Eq. (43), this factor substantially changes the coefficients ri^'^, 

%=3 , ^2 = -2 ' ^^4 = 3^, "^0=3 ' '^2 = -g , ''4 = -Y^- (64) 
This leads to d = 6. 

Wo finally comment on why it is impossible that a purely longitudinal order parameter induces a transverse gap. 
Inserting the matrix A^k from Eq. (59) into the anomalous propagator from Eq. (25), and the result into the 
right-hand side of the gap equation (26), we realize that the resulting Dirac structure still commutes with 75 and thus 
preserves the chirality. This is the characteristic feature of a longitudinal gap. Therefore, the ansatz (59) does not 
induce a transverse gap on the right-hand side of the gap equation. 
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F. The transverse CSL phcise 



For transverse gaps, 

Mk = J-7±(k). (65) 
The condition (17) is fulfilled because 7j_(k) commutes with the energy projector A^. For the matrix Lk we obtain 

(Lk)^^ = 2^'F<5„;3-7L^(k)7i^^(k) • (66) 
The eigenvalues of this matrix are (cf. Appendix A) 

Ai = 2 (8-fold) , A2 = (4-fold) . (67) 

The projectors 7-"^'^ are given by 



(68) 



Although J • 7_L(k) 7-"^ 7^ for both s = 1 and s = 2, the final result for T^'ff j'^(k, q) is nevertheless zero. To see this, 
however, one has to explicitly perform the trace in Eq. (28). For s = 1, we obtain 



Tm'^i^<l) = \ (l + ee't-q)' , 



(69a) 
(69b) 



The constants a^. are the same as in the 2SC and longitudinal CSL phases, see Eq. (43). The coefficients 773^ are 



%=3 



% = -3 



^^^=24 



%=3 



^^=24 



(70) 



This gives d = 9/2. 

For the same reasons as explained at the end of the last subsection, it is impossible to induce a longitudinal gap 
with the matrix A4k of Eq. (65) on the right-hand side of the gap equation. 



G. The polar phase 

In contrast to the CSL phase, in the polar phase the vector <p'^{K) in Eq. (51) does not couple to color space. 
Instead, it simply points into a fixed spatial direction, which we choose to be the ^-axis. Consequently, the matrix 
Mk in Eq. (16) reads 



fc^+7l(k) . (71) 



As in the 2SC case, the condensate is aligned with the (anti) blue direction in color space. Thus, condensation 
spontaneously breaks the color SU{3)c and spatial SO{3) symmetries to SU {2)c and SO{2), respectively. 

Due to the identity (k^ + j^){k^ — 7^) = 1, the Dirac structure of the matrix Lk is trivial, and it looks rather 
similar as in the 2SC case, Eq. (39), 

(ik)^^ = iJlT' = (S'' - S^') Sa0 . (72) 

This similarity is also apparent in the eigenvalues of Lk (cf. Appendix A), 

Ai = 1 (8-fold) , A2 = (4-fold) , (73) 

where the degeneracy refers to the combined color and Dirac spaces. The projectors are the same as in Eq. (41). For 
this reason, we again immediately conclude that 7^o^j'^(k, q) = 0. For s = 1 we obtain 
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To-''i(k,q) = i {(l + ee'k-q) 



Tr'''(k,q) = 3 hk^r (l-ee'k-q) 



1 + (1 + ee') y r] - (efc' + e'ff} 
{ek-e'qf 



1 



X { (l + ee' k • q) [l + (1 + ee') F ^ ] - (efc^ + e'r )^ }) 



(74a) 



(74b) 



Prom this it is obvious that oi = 1 and 02 = as in the 2SC case, cf. Eq. (43). 

There is, however, a marked difference between the expressions (74) and the corresponding ones for all previously 
discussed cases. In contrast to the other cases, there arc two independent, fixed spatial directions, that of the order 
parameter and that of the vector k. Since we already aligned the order parameter with the ^;-direction, we are no 
longer free to choose k = (0, 0, k) for the d^q-integration. Without loss of generality, however, we may assume k to 
lie in the x^-plane, i.e., k = fc(sint?, O,cos^?), where '& is the angle between the order parameter and k. In spherical 
coordinates for the d'^q-iutegration, the azimuthal angle 9 is no longer identical with the angle between k and q. Or 
in other words, k • q no longer depends solely on 9, but also on the polar angle ip. This has the consequence that also 
the modulus of the gluon 3-momentum p depends on ip. Since p enters the gluon spectral densities in a complicated 
fashion, it appears impossible to perform the (^-integration analytically in this way. 

The solution is to rotate the coordinate frame for the (i'^q-integration by the angle ■& around the y-axis, such that 
the rotated 2;-direction aligns with k. The quantities k • q, , and appearing in Eqs. (74) are expressed in terms 
of the new spherical coordinates (g, 9' , (p') and the rotation angle 1? as follows: 



k • q = cos 9' 



= cos 9' cos 1? — sin 9' sin 1? cos <p' 



k^ = cos 



(75) 



In the new coordinates the angle between k and q is identical with the azimuthal angle 0', and thus p becomes 
independent of (p' . Still, the v5'-integral is not trivial because of the potential p' dependence of the gap function. At 
this point we can only proceed by assuming the gap function to be independent of p' . With this assumption, the 
(^'-integration becomes elementary, and we are finally able to read off the coefficients f]2m' which now depend on 1?, 



%=3 



V2 



2 + cos^ •& 



1 + cos^ 'd 
24 



%=3 



2 - cos^ 'd 



1 -Scos^'d 
24 



(76) 



From this and Eq. (37) we compute d = 3(3 + cos^ i?)/2. 

Let us now comment on our assumption that the gap function is independent of (/?'. As mentioned in the introduction 
and as will be shown in the next section, the value of the gap function at the Fermi surface, (j)o, is proportional to 
exp(— d), cf. also Table I. The angular dependence of d then implies a similar dependence of the gap itself. If k 
points in the same direction as the order parameter, i9 = 0, we find d = 6, while for k being orthogonal to the order 
parameter, i} = ■jt/2, one obtains d = 9/2. In the first case, the gap is longitudinal in the sense introduced in Sec. II D, 
while in the second it is transverse. These two cases have also been discussed in Refs. [7,10], with the same results 
for the constant d. Our results surpass the previous ones in that they interpolate between these two limiting cases. 

However, the angular dependence of (f>o causes the following problem. The gap function (j){ek,i,^) is proportional 
to cf. the next section, and thus also depends on i?. Under the (i"^q-integral on the right-hand side of the gap 
equation, this dependence translates into a p' dependence of (j){eq^i, q). Our previous assumption, which was necessary 
in order to perform the (^'-integral, precisely neglected this dependence. Therefore, this approximation is in principle 
inconsistent. Nevertheless, the agreement of our results with the ones of Refs. [7,10] suggest that the (^'-dependence 
of the gap function could be a sub-subleading effect. 



III. SOLUTION OF THE GAP EQUATION 

In this section wc solve the gap equation (4). Let us first distinguish between the cases where oi = 1, 0,2 = 0, and 
where both ai and 02 are nonzero. The former are the 2SC phase. Sec. II B, the longitudinal and transverse CSL 
phases. Sees. HE and II F, and the polar phase, Sec. II G. The latter are the CFL phase. Sec. II C, and the CSL phase 
with both longitudinal and transverse gaps. Sec. II D. 

In the former cases, there is only one gapped quasiparticle excitation and the solution of the gap equation (4) is 
well-known. It was discussed in detail in Ref. [8]. In the 2SC phase, the longitudinal CSL phase, and the polar phase, 
all one has to do is replace the constant h in the calculation of Ref. [8] by the constant h = &exp(— d), cf. Eq. (6). 
The result for the value of the gap function at the Fermi surface is Eq. (7), but without the factor (A^^ A2^)~^/^. 
However, one immediately reads off Table I that in the respective cases this factor trivially equals one. In the 2SC 
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phase, d = 0, and consequently b = b, such that the result coincides with Eq. (1). In the other phases, where d > 0, 
the gap is reduced as compared to the 2SC phase by a factor exp(— rf), cf. Table I. 

There is a slight subtlety when solving the gap equation in the transverse CSL phase. The value of the nonvanishing 
eigenvalue is Ai = 2, not 1. One has to multiply both sides of Eq. (4) with ^/Xl in order to obtain a gap equation 
for which the solution of Ref. [8] applies. This rescaling is appropriate, as in this case the gap in the quasiparticle 
excitation spectrum is indeed 2^/Xl(j)o, and not simply 2(f>o- The factor (A"^ X2^)~^/^ in Eq. (7) precisely accounts 
for this rescaling of the gap function, such that this equation is also valid in the transverse CSL phase. 

In the CFL phase and the CSL phase with both longitudinal and transverse gaps, there are two gapped quasiparticle 
excitations, which renders the solution of Eq. (4) somewhat more complicated. A priori, one has to solve two gap 
equations, one for each quasiparticle mass shell, ko = Ch^i and kg = efc_2- Therefore, as a function of momentum k, 
there are in principle two different gap functions, (t)r{k) = <p{ek,r,k), r = 1,2. 

In order to proceed with the solution, to subleading order we may approximate the logarithm in Eq. (4) in a way 
first proposed by Son [6], 




c± e{eg,s - ek,r) In ( ) + e{€k,r - eq,s) In ( — ) ■ 

^o,s J \ ^k,r , 



With this approximation and the new variables 



= g In 



k-ii + tk 

to subleading order the gap equation (4) transforms into [8] 



2/s = In 



26/i 



q- fJ, + Cq^i 



\ f^^ civ )" f^^ 

) = ^aAxr dys{l-2g y^) tanh — ^ (f>{ys) + dysys{l-2g y^) tanh 



2T 



Here, we denoted the value of Xs at the Fermi surface, i.e, iov k = fi and ek,s = ^^,8, by 

2&/i 



X* = 5 In 



■S 90,8 



(77) 



(78) 



(79) 



(80) 



where 0o,s = <Pi.^l) is the value of the function (p{xg) at the Fermi surface. The single point k = jj in momentum 
space thus corresponds to two different points xl, X2, xl ^ X2 , in the new variables Xg. Since we expect (f)o,s to be 
~ exp(— 1/^), a;* is a constant of order one. Furthermore we defined 



xo=g\n[ — 



(81) 



This constant is parametrically of order 0{g). To subleading order, the relation between the new variable yg and the 
excitation energy is given by [5], 



(82) 



e(ys) = 6/x exp — - 



A consequence of the transformation of variables (78) and of neglecting sub- subleading corrections is that the two 
equations (79) for r = 1 and r = 2 become identical. The only difference is the notation for the argument of the 
function 0, which in both cases we may simply call x. Therefore, instead of two separate equations, we only have to 
consider a single equation which determines the function (^{x). Moreover, y^ is merely an integration variable, and 
we may set = y in the following. 

With Eq. (5), we rewrite Eq. (79) in the form 



J X 



dy {1 — 2gy) tanh 



2T 



■f. 



a\x I (1 — 2^7/) tanh 



0(y) 
<y) 



dyy{l-'2.gy) tanh 



<y) 



2T 



2T 



m ■ 



(83) 
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One can also write this equation in a form where X2 is replaced by and ai by a2, respectively. Equation (83) 
is an integral equation for the function 4>{x), which is solved in the standard manner by converting it into a set of 
differential equations [6], 



dx 
dx"^ 



= I dy{l -2gy) tanh 

■e{x) 



2T 



dy {1 — 2 gy) tanh 



<y) 



2T 



= —(1 — 2gx) tanh 



2T 



(x). 



(84a) 
(84b) 



We now solve the second-order differential equation (84b) at zero temperature, T = 0. One immediately observes that 
this equation is identical to Eq. (22c) of Ref. [8] , and its solution proceeds along the same lines as outlined there. The 
only difference compared to the previous calculation are the extra terms ai in Eqs. (83) and (84a). To subleading 
order, we expect 4'o,i/^o,2 — 1 (we show below that this assumption is consistent with our final result), such that the 
difference 



x^-xl=gln 



'0.1 



(85) 



is of order 0{g). Consequently, the extra terms ~ ai are of subleading order, O{g(j)o), and we may approximate 



dy{l-2gy) (p{y) ~ {x^ - xl) <j>o^: 



(86) 



Since wo always ordered the eigenvalues siicli that Ai > A2, cf. Table I, X2 — x^ > 0. 

The subleading correction (86) qualitatively changes the behavior of the gap function (j){x) near the Fermi surface. 
In the absence of the term ai in Eq. (84a), the derivative of the gap function vanishes for x = X2, and the gap 
function assumes its maximum at this point [8]. The subleading correction (86) induced by the two-gap structure 
in the CFL and CSL phases causes the derivative (84a) of the function (l){x) to be negative at the Fermi surface. 
Consequently, since we still expect (f){x) to rapidly vanish away from the Fermi surface, this function assumes its 
maximum not right at the Fermi surface, but at a point Xmax which is close, but not identical to X2- We shall see 
that X2 - a;,„ax ^ 0{g). 

The subleading correction (86) modifies the solution of the differential equation (84b) from the one given in Ref. 
[8]. Again, we fix the two unknown constants in the general solution of the second-order differential equation (84b) by 
matching the solution and its derivative to the right-hand sides of Eqs. (83) and (84a) at the point x = X2- Introducing 
the variables z = -{2g)~'^/^ (1 - 2gx) and z* = -{2g)~'^/^ (1 - 2^0:2), the solution reads 



'0,2 



\ M{\z\) smMz*\)-6{\z\)] 
1m(|^*|) sin [<^(|^*|)- ^(1^*1)] 



+ ai {x*2-xX){2g)-^'^ 



M{\z\) sm[e{\z*\)-e{\z\)] 

N{\z*\) sin [<^(|^*|)- ^(1^*1)] 



(87) 



where the functions M(|z|), A^(|z|), <^(|2;|), and 9{\z\) are related to the Airy functions Ai(z), Bi(z) and their derivatives 
in the standard way [21]. The derivative d(l){z)/dz can be obtained from Eq. (87) simply by replacing M(|2|) and 
0{\z\) by Ai'dzj) and <^(|2|), respectively. The difference to the solution for a single gapped quasiparticle excitation, 
cf. Eq. (27) of Ref. [8], is the term proportional to ai. 

Finally, we have to determine the value of (/'o,2- To this end, we rewrite Eq. (83) at the point a; = in the form 



zo + {2g)-y^] ^izo)=<l>{zo 



dz 



(88) 



where zq = —{2g)~^^^ (1 — 2gxo). Remarkably, this equation holds in this form also in the case of a single gapped 
quasiparticle excitation, cf. Eq. (29) of Ref. [8]. In weak coupling, the dependence on the variable zq is spurious. 

Inserting the solution (87) and its derivative for z = zq and expanding M(|2;o|), A^'d^oDi V'd^^oD) and ^(|-2o|) to order 
0{g) as demonstrated in Ref. [8], one derives the condition 



+ ai {X2 - xl) . 



(89) 



The second term is the 0{g) correction originating from the quark self-energy. It leads to the constant 6q in Eq. (1) 
and was first derived in Refs. [7,8]. The last term ai is the correction arising from the two-gap structure in the 
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CFL and CSL phases to the resuh (33) of Ref. [8]. Because of Eq. (85), this correction is also of order 0(g). Using 
the definition (80) of X2, as well as the condition (5), we conclude that the expression for (/)o,2 is identical to the one 
for 00 in Eq. (7). This is the value of the gap function at the Fermi surface, fc = /i, or x = Xg, for the quasiparticle 
excitation branch €k,2- The additional suppression factor compared to the 2SC gap (pQ^''' of Eq. (1), which originates 

from the two-gap structure, is (A^^ Aj^)"^^^. For the CFL phase, we obtain the value 2"^/-^, while for the CSL phase, 
we have 2"^/^. 

We can also compute the gap fimction at the Fermi surface for the first excitation branch ek,i, i.e., at a; = a;|. The 
difference 0o,2 — (po,i can be obtained from Eq. (83) as 

-^0,2 - <^o,i = f ' dy[y-xl-ai {x*^ - x^)] {l-2gy) cl){y) . (90) 

An upper bound for the term in brackets is given by setting y — where it assumes the value 02 (xj — x\) on account 
of Eq. (5). Pulling this factor out of the integral, the latter can be estimated with Eq. (86). This proves that the 
difference ^0,2 — </'o,i is only of order O{g^4'o), which shows that our above assumption (/)o,i/0o,2 — 1 is consistent up 
to subleading order. To this order, we may therefore set 0o,i = </'o,2 = (j^o- 

We now determine the value of Xmax, where the gap function assumes its maximum, by setting the left-hand side 
of Eq. (84a) equal to zero. This leads to the condition 

r dy{l-2-gy) cj>{y) = ai H dy{l-2gy) cj>{y) . (91) 

To order O{g4>o), one may easily solve this equation for Xmax, with the result 

a;max = Xs - «i I In (^^^ ) (92) 

i.e., Xmax is indeed smaller than Xj by a term of order 0{g), as claimed above. Obviously, since ai < 1, from Eq. 
(85) we derive the inequality Xj < Xmax < •^21 i-^-i the gap hmction assumes its maximum between the values x^ and 
X2. The value of the gap function at Xmax can be estimated via a calculation similar to the one for the difference 
^0,2 — ^0,1 above. The result is ^moa; — ^o[l + 0{g'^)]. This means that the gap function is fairly flat over a region 
of size 0{g) (in the variable x) in the vicinity of the Fermi surface. 



IV. TRANSITION TEMPERATURE 



In this section we compute the transition temperature Tc where the color-superconducting condensate melts. In 
the 2SC phase, in the CSL phase with longitudinal gap, and in the polar phase, the calculation of Ref. [8] applies, 
and we obtain the BCS result Tc/(f)o = e'^/n. In these cases, (A^^Aa')^/^ = 1, cf. Table I, such that Eq. (9) is valid. 

In the transverse CSL phase, we also obtain the BCS result for the relationship between Tc and the zero-temperature 
gap after a rescaling of the gap function by a factor \/Xi, cf. the discussion in Sec. III. Since in this case (A^^Aa^)^/^ = 
\/Ai, Eq. (9) also applies. 

In the CFL phase and the CSL phase with both longitudinal and transverse gaps, wc have to compute Tc explicitly. 
The calculation follows the line of arguments presented in Ref. [8], taking into account the additional term ~ oi in 
Eq. (83). As in Refs. [5,8] we assume that, to leading order, the effect of temperature is a change of the magnitude 
of the gap, but not of the shape of the gap function, 



cj>{x,T)c^cj){T) 



'/'(^,0) 



(93) 



where (f){T) = 4'{x2, T) is the value of the gap at the Fermi surface at temperature T, 0(x, 0) is the zero-temperature 
gap function <p{x) computed in the last section, cf. Eq. (87), and = (^0,2 = 0(a;2,O). With this assumption, Eq. 
(83) reads at the Fermi surface 



1=1 dyy{l - 2gy) tanh 

J XO 



2T 



-^(y.o) 



dyy{l - 2gy) tanh 



2T 



-ai X2 / dy{l-2gy) tanh 
; J1+J2+J3 , 



2T 



(94) 
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where we divided the second integral in Eq. (83) into two integrals: Xi which runs from xq to a;^, with x„ = 
X2 — g ln(2K), K ^ 1, and X2 which runs from x,^ to X2 [5]. We now compute the integrals Ii through I3 separately 
to subleading accuracy, i.e., to order 0{g). 

In the first integral Ji, which runs over a region far from the Fermi surface, e(y) T, and we may approximate 
the tanh by 1. This integral can be formally solved by integration by parts using the differential equation (84b), 



00 



0(a;K,O) - ^(a;«,0) 



(95) 



where we exploited the condition (88). Expanding the functions on the right-hand side around X2 we obtain to 
subleading order 



Ji = 1 - - [g ln(2K) - ai {x*2 - x\)] 



(96) 



This estimate is similar to the one made in Eq. (36) of Ref. [8]. The main difference to that calculation is the term 
~ tti which appears because the first derivative of the gap function no longer vanishes at the Fermi surface, cf. the 
discussion in the previous section. 

In the second integral I2, which only contributes to order 0{g) to the right-hand side of Eq. (94), to subleading 
order we may set (j){y, 0)/(/)o — 1 and y ~ a;^ — 7r/2. Reverting the transformation of variables (78) we obtain 



^2 = TTff 



^"^o d{q - n) 



tanh 



\2T ) 



(97) 



The last integral in Eq. (94), I3, also contributes a term of order 0{g), and may thus be approximated by an 
argument similar to that leading to Eq. (86), 



X3 = a\ X2{x2 — x\) tanh 



4>{T) 
2T 



(98) 



At the critical temperature T^., where </>(Tc) 
becomes 



0, this term vanishes. Putting everything together, at T = Eq. (94) 



9 1 d{q 



1^) 



1 



tanhf^ 



-ai {x*2 - xl) 



(99) 



where the term ln(2K) in Eq. (96) was expressed in terms of an integral according to Eq. (96) of Ref. [5]. In the 
integral on the left-hand side, we may send k ^ 00 [5]. This allows us to perform it analytically, which yields the 
result \n.[e^ ^/\2(f)Q / {'kTc)]^ where 7 ^ 0.577 is the Euler-Mascheroni constant. If the right-hand side of Eq. (99) were 
zero, for A2 = 1 this would then lead to the BCS relation Tc/4>o = e"' /tt. However, using Eq. (85) we now obtain 
Eq. (9). The last factor on the right-hand side of this equation is exactly the inverse of the additional factor in Eq. 
(7). This factor violates the BCS relation Tc/(j)o = c'^/tt in the CFL and CSL cases. In the first case, the transition 
temperature is by a factor 2^/^ larger than one would expect from BCS theory and in the second case it is larger by 
a factor 2^/"^. However, in units of energy, this factor just cancels the one from 0o in Eq. (7). 
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APPENDIX A: COMPUTING EIGENVALUES 



The eigenvalues A^. of Lk follow from the roots of 

det (A 1 - Lk) = . 
The left-hand side of this equation can be rewritten in the form 



(AI) 
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det ( A 1 - Lk) = exp {Tr [In ( A 1 - Lk)] } . ( A2) 

The logarithm of the matrix A 1 — Lk is formally defined in terms of a power series, 



Tr [In (A 1 - Lk)] = In ATr 1 + Tr 



InM-^J =lnATrl-^-A-"TrXi^ . (A3) 

^ ' -' n=l 



In order to proceed, one needs to know the trace of the nth power of the matrix Lk- In the 2SC phase, the CSL 
phase with longitudinal gap, and in the polar phase, this is particularly simple, since Lk is a projector, cf. Eqs. (39), 
(60), and (72), hence = L^. Counting color and flavor degrees of freedom in the 2SC phase, and color and Dirac 
degrees of freedom in the longitudinal CSL and polar phases, the trace of ik is 4 in the former and 8 in the latter 
case, respectively. Therefore, we obtain for the 2SC phase 

2SC phase: det (A 1 - Lu) = A^ (A - 1)^ = . (A4) 

This yields the eigenvalues given in Eq. (40). For the longitudinal CSL and polar phases we analogously compute 

long. CSL and polar phases: det (A 1 - Lk) = A"* (A - 1)^ = , (A5) 

which leads to the eigenvalues of Eqs. (61) and (73). 

The next simple case is the CSL phase with transverse gap. In this case, Lk is not a projector, but since 7j_-7_l = —2, 
it is still idempotent up to a factor, = 2Lk- Because of LJJ = 2"~^ Lk and TrLk = 16 we then obtain 

trans. CSL phase: det (A 1 - Lk) = A^ (A - 2)^ = , (A6) 

from which we read off the eigenvalues of Eq. (67) . 

The CFL phase and the CSL phase (with both longitudinal and transverse gaps) are the only cases where the 
calculation of is slightly more involved. First, one proves the identity L^ = 5 Lk — 4 1, which is valid in both cases. 
Repeated application of this relation allows to reduce an arbitrary number of powers of Lk to a single power, plus a 
term proportional to the unit matrix, 

LjJ = a„ Lk + 6„ 1 . (A7) 
Multiplying both sides of this equation by Lk, one derives the recursion relation 

a„+i = 5a„ - 4a„_i (AS) 

for the coefficients a„, and the identity 

b„+i = -4 a„ (A9) 

for the coefficients 6„. The recursion relation (A8) can be solved with the Ansatz a„ = p", which yields a quadratic 
equation for p with the solutions pi = 4 and P2 = 1. The general solution of the recursion relation is then a„ = 
api + (3p2 = a4" + /?. The coefficients a and /3 can be determined from ai = 1 and 02 = 5, such that 

4" _ 1 4" — 4 
an = , bn = — . (AlO) 

In the CFL phase, TrLk = 12 and Tr 1 = 9, while in the CSL phase, TrLk = 24 and Tr 1 = 12. Consequently, in the 
CFL phase 

CFL phase: det (A 1 - Lk) = (A - 4) (A - 1)^ = , (All) 
which leads to Eq. (47), while in the CSL phase 

CSL phase: det (A 1 - Lk) = (A - 4)^ (A - 1)*^ = , (A12) 

which yields Eq. (54). 
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APPENDIX B: INTEGRATION OVER GLUON MOMENTUM 



In this appendix we compute the integrals over gluon 3-monientuni p to subleading order in the gap equation. We 
shall see that to this order it is consistent to put k — q — fi. 

After replacing k • q = (fc^ + — p^)/ (2kq), the coefficients Tii2m{^^' ^ ^> can be read oft from Eqs. (42), (49), (56), 
(63), (69), and (74). One first observes that for all cases considered here, ??2m(s6'> ^) = foi' w > 3. Next, one also 
realizes that r]^_2 = 0, since there is no term in Tqq proportional to Consequently, we have to compute the 

integrals 



k+q 2 



for the contribution of static electric gluons to the gap equation. 



k+q 



M 



dpp y- 
\kq 



m = 0,1,2 



m = -1,0, 1,2 



(Bla) 



(Bib) 



for the contribution of non-static magnetic gluons, and 



j-t,2 
'-2m 



J\k-q\ 



dpp 



P 



p6 + M'^ul 



m 



-1,0,1,2, 



(Blc) 



with uj± 
(Bib) is 



g ± e% ^, for the contribution of almost static magnetic gluons. The result for the integrals (Bla) and 



{k + qf + 3m 



(fc - qf + 3m2 



m 



kq 



k^ + q- - Wml 



kq 



rt,l 



rt,l 



-^4 



kq 
In 



kq 



{k + qf 
{k + qf ^ 
M2 _ " 

(fc + qf - 

kq 

{k + qY 




2{kqy 



(B2a) 
(B2b) 

(B2c) 

(B2d) 
(B2e) 
(B2f) 
(B2g) 



The approximate equalities on the right-hand sides hold to subleading order in the gap equation. One obtains 
them employing two approximations. First, terms proportional to at least one power of or carry at least two 
additional powers of g, which renders them sub-subleading and thus negligible to the order we are computing. Second, 
one utilizes the fact that the q integration in the gap equation is over a region of size 25 around the Fermi surface, 
where 5 ^ nig. To subleading order it is thus accurate to put k = q = ^ (see discussion in Sec. II A). This then yields 
the right-hand sides of Eqs. (B2). 

Note that there is a term ~ /x^/M^ ~ in Eq. (B2d). This term is parametrically the largest and could in 
principle give the dominant contribution to the gap equation. However, in all cases considered here, it turns out that 
the coefficient ryij is proportional to at least one power of {k — q)^. Performing also the q integration in the gap 
equation, one then has terms of the form 



Jo 



d{q -fl) {k- qf 



M2 



4>{<^q,Q) 



(B3) 
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where for the purpose of power counting we have neglected the q dependence of the gap function, (f){eq, q) ^ (j}^, and 
we have evaluated the integral on the left-hand side for fc = /i. As long as (5 ~ Wg ~ M, the leading term in Eq. (B3) 
is ^ g^i/'o, and thus it is only of sub-sublcading order in the gap equation. It is obvious that the constant term —1/4 
in Eq. (B2d) is parametrically even smaller. The contribution to the term ~ rf_2 from non-static magnetic gluons is 
therefore negligible to subleading order. 

Finally, also the integrals can be computed analytically [22]. Defining a = (M^w^)^/^, the result is 



-7-4,2 
-^-2 



12a \ 



In 



— ax + 



— 2\/3arctg 



2x — a 



{k-qY 



rt,2 



7-t,2 



-^4 



6 



M6 + a3 



{k - q)^ + 
M^-{k- qf 











6 \ 





2kq Ukq 

M4 - (fc - qy 



In 



{x + af 



x^ — ax + a^ 



2^3 arctg 



2x — a 
\/3a 



{k-qY 



0, 



A{kqY 



(kq) 



(B4a) 

(B4b) 
(B4c) 
(B4d) 



Here, we used the short notation {f{x)}^ = /(a) — f{b). In order to obtain the approximate equalities on the right- 



hand sides of Eqs. (B4b), (B4c), and (B4d), one employs the fact that typically (fc — q)"^ 
parametrically {k — qY <C a <C Af^. This immediately yields the right-hand side of Eq. (B4b) 



< M^, such that 
For Eqs. (B4c) and 



(B4d), we use this estimate in order to expand the logarithm occurring in Eqs. (B4a) and (B4c). One finds that 
the leading term is a/M'^. Similarly, one expands the inverse tangent occurring in these equations, which leads to 
terms which are even of order 0(1). Collecting all prcfactors, however, all terms in Eqs. (B4c) and (B4d) are then 
suppressed by at least one power of . These sub-subleading corrections arc negligible to the order we are computing. 

Somewhat more care is necessary in estimating the; tcirms in Eq. (B4a). Again, one may expand the logarithm and 
the inverse tangent. Together with the prefactor, this leads to a term ~ l/M"^ for the logarithm, and a term r-^ 1/a 
for the inverse tangent. The first term is harmless: together with the factor (fc — qY from 77^.2 it leads to an integral 
of the form (B3), which was already shown to give a sub-subleading contribution to the gap equation. The other term 
leads to the integral 



Jo 



d{q -11) (k- qf 



a 



H^q,q) 



where we used similar power-counting arguments as in Eq. (B3). The last integral is finite even for 



(B5) 



0, so that 



this contribution is then again ^ 5 (/>o and thus of sub-subleading order 



we can estimate it to be ^ S^^'''. For S 
in the gap equation. 

In conclusion, also the contribution of almost static magnetic gluons to the term ~ t]*_2 is of sub-subleading order 
and can be neglected. To subleading order, it is therefore consistent to put 77^3 = from the beginning, provided one 
chooses 6 rUg. 



APPENDIX C: DIRAC STRUCTURE OF THE SPIN-1 GAP MATRIX 

Since the gap matrix T,2i{K) is a complex 4x4 matrix in Dirac space, it can be written as a linear combination of 
sixteen basis matrices. The gap matrix is a scalar in momentum space, which reduces this number to eight [12]. We 
choose this basis set to be 

M = (Ml, ...,Ms) = (1,7°,7 • k,7 • k7°,7 • k7V,7V,7 • k7',7') • (CI) 

In the spin-1 case the order parameter has to be a three-dimensional vector. But E21 has no vector structure. Thus 
the order parameter has to be contracted with other vectors. The only available vectors are k and 7, and each 
contraction can still multiply any element of M. Thus, we can write the gap matrix in terms of 16 3-vector order 
parameters lp^, i = 1, . . . , 16, as 

8 16 

S2i(if) = J2'^i(^) ■ + E V'iW • 7Mi_8 . (C2) 
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Decomposing the ith vector order parameter ip^ into a longitudinal and a transverse part with respect to k, ^p^ = 
+ (fil, where y;- = • (1 — kk), the gap matrix becomes 



8 16 



Mi. 



(C3) 



i=l 



Note that whenever a basis matrix Mj is multiplied (from the left) by 7-k one obtains, up to a minus sign, another basis 
matrix Mj. Thus we can rearrange this expression such that there are only eight different longitudinal coefhcients. 
The eight longitudinal and sixteen (independent) transverse coefficients can be combined into eight new 3- vector order 
parameters 



S2i(if) = ^0,(if)- [k + 7±(k) 



Mi 



(C4) 



where 7j_(k) = 7 • (1 — kk). Now one can perform a basis transformation and write the gap matrix in terms of the 
projectors for energy, chirality and helicity. This is completely analogous to the spin-0 case [12]. In the ultrarelativistic 
limit two out of the three projectors arc sufficient, for instance those for energy, A|, and chirality, Vh = (1 + hj^)/2, 
h = ± for right- or left-handed quarks, respectively. Then we are left with four 3-vector order parameters in terms of 
which the gap matrix reads 



(C5) 



e,h 



For condensation in the even-parity channel, only two 3-vector order parameters are independent, because (f)^ — <f}'^ = 
(jf [5]. The sum over chirality projections can be immediately performed to give 



E2i(i^) = 5]</)^(/^)- [k + 7^(k) 



A 



k • 



(C6) 



This is Eq. (51). 



[1] D. Bailin and A. Love, Phys. Rep. 107, 325 (1984). 

[2] M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. B 422, 247 (1998); 

R. Rapp, T. Schiifer, E.V. Shuryak, and M. Velkovsky, Phys. Rev. Lett. 81, 53 (1998). 
[3] D.J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973); H.D. Politzer, ibid. 1346. 

[4] T. Schafer and F. Wilczek, Phys. Rev. D 60, 114003 (1999); D.K. Hong, V.A. Miransky, LA. Shovkovy, and L.C.R. 

Wijewardhana, Phys. Rev. D 61, 056001 (2000); Erratum ihid. D 62, 059903 (2000). 
[5] R.D. Pisarski and D.H. Rischke, Phys. Rev. D 61, 074017 (2000). 
[6] D.T. Son, Phys. Rev. D 59, 094019 (1999). 

[7] W.E. Brown, J.T. Liu, and H.-C. Ren, Phys. Rev. D 61, 114012 (2000); 62, 054013 (2000); 62, 054016 (2000). 

[8] Q. Wang and D.H. Rischke, Phys. Rev. D 65, 054005 (2002). 

[9] T. Schafer, Nucl. Phys. B575, 269 (2000). 
[10] T. Schafer, Phys. Rev. D 62, 094007 (2000). 
[11] R.D. Pisarski and D.H. Rischke, Nucl. Phys. A702, 177 (2002). 
[12] R.D. Pisarski and D.H. Rischke, Phys. Rev. D 60, 094013 (1999). 
[13] M. Alford, K. Rajagopal, and F. Wilczek, Nucl. Phys. B537, 443 (1999). 
[14] M. Iwasaki and T. Iwado, Phys. Lett. B 350, 163 (1995). 

[15] D. VoUhardt and P. Wolfle, The Superfluid Phases of Helium 3 (Taylor & Francis, London, 1990). 

[16] J. Bardeen, L.N. Cooper, and J.R. Schrieffer, Phys. Rev. 108, 1175 (1957); M. Tinkham, Superconductivity (Gordon and 

Breach, Science Publishers, New York, 1965). 
[17] C. Manuel, Phys. Rev. D 62, 114008 (2000). 
[18] D.H. Rischke, Phys. Rev. D 64, 094003 (2001). 
[19] K. Zarembo, Phys. Rev. D 62, 054003 (2000). 



21 



[20] R.D. Pisarski and D.H. Rischke, Proc. of A Symposium in Memory of Judah M. Eisenberg: Nuclear Matter, Hot and Cold, 
Tel Aviv, Israel, April 14 - 16, 1999, edited by J. Alstcr and D. Ashery (Tel Aviv University Press, Tel Aviv, 2000), p. 
101; e-print archive nucl-tli/9907094. 

[21] Handbook of Mathematical Functions, edited by M. Abramowitz and LA. Stegun (Dover, New York, 1965). 

[22] I.S. Gradsteyn and I.M. Ryzhik, Table of integrals, series and products, (Academic Press, San Diego, 1990), Eqs. (2.126.2)- 
(2.126.5). 



22 



